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Beschreibung 




Die Erfindung betriff t ein Veifahren zur Abbildung eines techniscben Systems durch tinea Rechnen Es ist wichtig, daB 
vorgebbar technische Systeme durch eine Modelliemng geeignet genau beschrieben werden. So ist es von groBer Bedeu- 
5 tung, ein komplexes technisches System, beispielsweise einen chemise hen ProzeB, der von mebreren EingabegroBen und 
mefareren AusgabegrSBen abhangt, geeignet genau zu modellieren (abzubilden), wobei ein Modellierungsfehler, also ein 
Uaterschied zwischen der Modelliemng und dem zugrundeliegenden realen System, minimiert werdeo solL 

Aus [1] ist bekannt, ein Differenti algleichungs system der Variationsgleichungeo anhand von Sensi ti vitatsmatrizen zu 
loseru Dieses Vorgehen weist deD Nachteil auf, daB die Berechnung ganzer Sensitivititsmatrizen viel Rechenzeit erfor- 
10 dert und daS damit ein auf dieser Methode beruhendes Verfahren deutliche GeschwindigkeitseinbuBen aufweist 

Aus [2] und [3] sind Verfahren der Lagrange-Multiplikatoxen bekannt Hieibei weiden flber eine LSsung der soge- 
nannten adjungierten Gleichungen die Lagiange-Multiplikatoren berechnet Jeder Mefipunkt liefert einen eigenen Bei- 
trag zu den Lagrange-MuMplikatoien, somit sind die adjungiexten Gleichungen projeweiligen MeBpunkt zu losen. Solch 
ein Verfahren ist sowohl auBerst zeitaufwendig, als auch stellt es hohe Anforderungen an den Speicherplatzbedarf auf 
15 dem Rechnex 

Aus [4] ist bekannt, Optimierungsalgorithmen fur groBe Systeme zu untersuchen. Dabei werden mit duTerentiellen 
und algebraischen Gleichungen beschriebene chemische Prozesse mit orthoganaler Kollokation diskretisiert und mit ei- 
ner sequentiell quadratischen Programmierung optimiert 

Die Aufgabe des erfindungsgemaBen Verfahrens besteht darin, einen Modellierungsfehler bei der Abbildung eines rea- 
20 len Systems zu minimieren. 

Diese Aufgabe wild gemaB den Merkmalen des Patentanspruchs 1 und des Patentanspruchs 2 gelttst 

Hierbei sei angemerkt, daB nachfolgend fettgedruckte Symbole in den Fbrmeldarstellungen \fektoren kennzeichnen. 

Der Modellfehler, der bei der Abbildung eines realen Systems in ein Modell entsteht wird dadurch beschrieben, daB 
eine B e wertungsfunktion definiert wild. Diese Bewertungsfunktion stellt die Abweichung der im realen System gemes- 
25 senen Ausgangswexte von den im Modell resultierenden Ausgangswerten dan Eine Optimierung der Beweitungsfunk- 
tion wild dadurch vorgenommen, daB auf dem Rechner Gradienten, die lokale Ableitungen der Bewertungsfunktion an 
vorgebbaren Werten von Systemparametem und einem Anfangszustandsvektor darstellen, ermittelt werden. Mit Hilfe 
der Gractientenfiinktion wild die Bewertungsfunktion minimiert und das System wird durch die Abbildung (das Modell) 
optimiert dargestellt 

30 Eine Weiteibildung des erfindungsgemaBen \ferfahrens besteht in der Verwendung einer quadratischen Abweichung 
als Bewertungsfunktion. 

Ferner ist eine Weiteibildung der Erfindung die Zerlegung eines zu betrachtenden Intervalls in mehrere Teilintervalle, 
auf die jeweiis das erfindungsgetnaBe Verfahren separat angewandt wild und beiechnete Tbilergebnisse zu dem Gesaml- 
ergebnis aufsummiert weiden. 
35 Die Erfindung wird anhand von AusfiJhrungsbeispielen, die in der Figur dargestellt sind, naher erlautert 

Es zeigt Fig. 1 ein Blockdiagramm, das Schritte des erfindungsgemaBen Verfahrens enthSlt 

In Fig, 1 weiden Schritte des erfindungsgemaBen Verfahrens dargestellt In video techniscbeo Anwendungen ist die 
Aufgabe gestellt, ein System (Ein-Ausgabe-System) der folgenden Form zu modellieren 



40 ^»<t)- fWtWt).w) 

(1), 



y(t) = m(x(t), w) 



wobei ein Modellierungspioblem in der Bestimmung der Systemparameter w besteht Gleichung (1) ist im Rahmen der 
Systemtheorie als Zustandsraumdarstellung eines nichtlinearen Systems bekannt Im Zustandsvektor x ist der interne Zu- 
stand des Systems dargestellt Eine zeilliche Entwicklung der internen Zustande wild durch ein dynamisches System (ein 

50 System von DifFerentialgleichungen) f beschrieben. Neben den konstanten Systemparametern w wild die Zeitentwick- 
lung ggf. noch durch externe GroBen, bier Eingabewerte u, die in der Regel zeitabhangig sind, beeinfluBt Ausgabewerte 
y des Systems, die das Systemveihalten nach auBen bin quantifizieren, hangen fiber eine statische Abbildung m in einer 
vorgebbaren Art und Weise vom Zustandsvektor x ab. Femer enthalt die statische Abbildung m im Normalfall Parameter; 
die ebenfalls in den Systemparametem w enthalten sind. 

55 Das Modellieiungsptoblem besteht in der Bestimmung der Systemparameter w des Systems, das durch Verbindung 
des dynanrischen Systems f mit der statischen Abbildung m beschrieben wird. Bekannt sind gemessene Ausgangsweite 
des (realen) System r^ zu vorgebbaren Zeitpunkten tv. Diese Zeitpunkte ^ sind monoton geoidnet und befinden sich in 
einem umfassenden Zeitintervall [to, T], inneihalb dessen das Systemverhalten untersucht wird. Die gemessenen Aus- 
gangsweite des System nv, die bei einem bestimmten Anfangszustandsvektor Xq die zeitliche Entwicklung der Ausgabe- 

60 werte der Modellierung y(t) in diskreter Form beschreiben, werden als eine MeBreihe bezeichnet. 

Bei Vorhandensein mehrerer MeBreihen sind alle im folgenden betrachteten Grofien fur jede MeBreihe einzeln zu be- 
stinimen und die Ergebnisse aufeusummieren. Nachfolgend wird also aus Griinden der Einfachheit und besseren Dar- 
stellbarkeit die Existenz nur einer MeBreihe angenommen. 

Besitzt das System Eingabewerte u(t), so sind diese ebenfalls als in ihrem Zeitvedauf vorgegeben angenommen. 

65 Das Modellierungsproblem ist gelost, wenn durch entspiechende Wahl der Systemparameter w in dem dynamischen 
System f und der statischen Abbildung m eine Obereinstimmung zwischen den gemessenen Ausgangswerten des Sy- 
stems n v und den zugehorigen Ausgabewerten der Modellierung yOJ erreicht wird. Die Ausgabewerte der Modellierung 
y(tv) stammen aus der Modellierung des (realen) Systems nach Gleichung (1) mit den berechneten Systemparametem w 
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im Zeidntervall [t^, T], au^^Bflfcran dem Anfangszustandsvektor Xq. Bei un^^^^wn Anfangszustand ist der An- 
fangszustandsvektor x<) eiir^B^^B'arameter des Modells und im Rahmen der^^PB&ung zusfitzlich zu den System- 
parametem w zu bestimmeiL 

Die techniscbe Losung des Modellieningsproblems wild dadurch erreicht, daB eine Optimierun gsaufgabe furmuliert 
wird. HIerzu wird fiber eine Bewertungsfunktioo die Abweichung zwischen gemessenem Ausgangswert des Systems Qy 
und Ausgabewerten der Modellierung y(tv) quantifizierL 

In der Praxis ist eine wichtige beispielhafte Form fiir die bewertete Abweichung gegeben durch eine quadratische Ab- 
weichung 

e(y(tv), n v ) = | (y(tv) - n v ) 2 (2). 

Allgemein setzt sich die BewertiingsfunktioD zusammen aus der Summe alls: bewerteten Abweichungen gemaB 

e = Z e(y(t v ), n v ) (3) - 



10 



15 



20 



die im Fall der quadrarischen Abweichung (siehe Gleichung (2)) auch als Fehlerquadratsumnte bezeichnet wild. 

Da die Bewertungsfunktion £ von den bereebneten Ausgabewerten der Modellierung y in Gleichung (1) abhangt, be- 
steht damit eine indirekte AbhSngigkeit zwischen der Bewertungsfunktion E und den Systemparametern w und ggf. von 
dern Anfangszustandsvekior x 0 . Eine Mimrnierung der Abweichungen zwischen den gemessenen Ausgangswerten des 
(realen) Systems und den beiechneten Ausgabewerten der Modellierung y^) entspricht einer Optimiening der Be- 25 
wertungsfunktion £. Das Modellierungsproblem ist damit Equivalent zu einer Optimierung der Bewertungsfunktion B 
durch entsprechende Variation der Systemparameter w und des Anfangszustandsveklors x©. Die Optimierung zur Be- 
stimmung von den Systemparametern w basiert auf einer Auswertung eines Gradienten 

— {Ableitung nach den Systemparametern) (gl) 30 

dw 

und bei unbekanntem Anfangszustandsvektor xq zusatziich auf einer Auswertung eines Gradienten 
rJE 35 

(Ableitung nach den Anf angs zu s t anden ) (g2) . 

dxo 

Das erfindungsgemaBe Verfahren eraioglicht die effiziente Bestimmung der beiden oben dargesteilten Gradienten (gl) 
und (g2), wobei pro MeBreihe nur eine Integration eines Systems von Differentialgleichungen anfallL 40 

Hierbei sei angemerkt, daB das Intervall [t^ T] alternativ in mehrere Teilintervalle zErf alien kann, wobei fiir jedes die- 
ser Ibilintervalle das erfindungsgemaBe Verfahren angewandt wird und die fur jedes lHlintervall ermitlelten GrdBen fur 
das Gesamtergebnis aufsummiert werderL 

Nacbfolgend werden die Schxitte der Erfindung zur Bestimmung der Gradienten dargestellL 

Im Schritt la wird das dynamische System 45 

A ^ = fWt) , ^ v ) (4) 

at 

vorwarts von einem Zeitpunkt to bis zu einem Zeitpunkt T unter Verwendung des Anfangszustandsveklors Xq, der be- SO 
kannt ist oder als zusatzlicher Parameter eingefflhrt wird, gelost Es ergibt sich ein Endzustand x(T). 1st der Endzustand 
x(T) anderweitig bekannt, so eatflllt Schritt la. 

Im Schritt lb wird ein System von Differentialgleichungen 

55 
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dt dw w 
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_ / 3f(x(t), u(t), w)V dS 
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d dB , . = ( dfjxjt), u(t), w)V 
dt dx 0 V J 



dE 
dXQ 



(t) 



(5) 



£a(t) o f(x(t), u(t), w) 

riickwarts von dem Zekpunkl T bis zu dem Zeitpunki to gel6st, wobei mit (.)' eine transponierte Matrix bezeichnet wild. 
Jedesmal bei Erreichen des MeBwertes iiy (Iq < < T) wild das LosungsvKfahren gestoppt (das gilt sinngemaB auch 
fur den Sonderfall ty = T) und nacb ciner diskontinuierlichen Anderung gem2B 



dx 0 v w dx 0 V VJ 



r 



dm 



dx 



V 



0e 



(6) 



wieder fortgesetzt bis zum nachsten MeBwert oder bis zum Erreichea des Anfangszeitpunkts (siehe Schritt lc). 

Fur den in Gleichung (2) genannten Spezialfall der Fefalerquadratsumme gilt anstelle von Gleichung (6) folgende An- 
derungsvorschrift. 



30 



35 



40 



dE / v dE / \ 

_ (ev)M __ (tv)+ 



f 

dm 




Jdx. 


«(tv); 



(m(x(t v ), w) - n v ) (6a) 



Als Anfangswertc einer Losung am Endpunkt T des Intervalls [to, T] sind dabei 

dB , v dE 

(T) : = 0 , 

dw V 7 dXQ 



(T):= 0, x(T) 



(7) 



zu vexwenden. 

In einem Schritt Id werden die Gradienten nach folgender Ybrschrift berechnet: 
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dm 
dw 



de 



dy(t v ) 



dE 



dE 



dac 0 dx 0 



(to) 



Fur den Spezialfall der Fehlerqu adratsu mme gilt dann: 

i 



- ?"?(*o) + Z 

dw dw ' 



dm, 



dw 



(m(x(t v ), w) - n v ) 



65 



dB 



dE 



(to) 



(8) . 



(8a) . 



dx 0 dbc 0 

Vorteilhaft an der Erfindung ist im Vergleich zur Methode der Varialionsgleichungen eine Performanz-Steigerung urn 
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eineo Faktox dim x. t^^^fefcv 

Ein weiterer Vorteil aer HB^^pbesteht darin, daS die Integration mit AusvS^P^pier Zwiscbcneigcbrdssc nur an 
den MeBpunkteo afoiderlichi^^Bt aber in einem zur numerischen Integialberechnung sehr fein zu wghlenden Raster. 

Weiterhin weist die Erfindung vorteilhaft die Bigenschaft auf, daB die Zahl der Integrationen eines DuTerentialglei- 
chungssystems mil der Anzahl der MeBreihen, nicht aber mit der Anzahl der MeBpunkte stedgt 5 

SchlieBlich ist das erfindungsgemaBe Yerfahren wesentlich genauer als eine Methode der finiten Differenzen. 
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Patentanspriiche 20 

1 . \ferfahren zur Abbildung eines technischen Systems durcb einen Rechner, 

a) bei dem eine Abweichung von gernessenen Ausgangswerten des technischen Systems und von Ausgangs- 
werten der Abbildung durch folgende Bewertungsfunktion definiert wild: 



E = Z K^v)^ n v) f 



25 



mit 30 



35 



^ x(t) = f(x(t), u(t), w) und 
y(t) = m(x(t),w), 

40 

wobei 

t eine Zeitvariable, 

y Ausgabewerte derModellierung, 

w Systemparameter, 

x einen zeitabhangigen Zustandsvektor; 45 

f0 eine zedtliche Entwicklung der internen Zustande der Abbildung des technischen Systems, 

u Eingabewerte, 

mO eine statische Abbildung, 

E die Bewermngsfunktion 

e eine bewertete Abweichung, 50 
v eine Zahl variable, 

Dy einen gernessenen Ausgangswert des Systems bezeichnen, 

b) bei dem durch den Rechner die Beweztungsfunktion dadurch optinnieit wird, daB Gradienten 

^ und 

dw dxo 

in folgenden Schritten bestimmt werden, wobei durch Xq ein Anfangszustandsvektor zum Zeitpunkt to darge- 
stellt wird: 60 
(1) das DirTerentialgleichungssystem 

A. x(t) = ^ ^ w) 

65 

wird von t = t$ bis t = T gelost, unter Verwendung des Anfangszustandsvektors wobei sich ein Endzu- 
stand x(T) ergibt, 
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Difiereatialgleichungen 



d dE 
dt dw 



(t) = - 



r d£(x(t), u(t), wfl dE (t) 



_d_ dE = _T af(x(t)/ u(t)/ w)^| dE x v 
dt dx 0 T cbc J dxo 



dt 



x(t) = f(x(t), u(t), w) 



wobei 0' die transponierte Matrix ist, wind von T riickwarts bis gelost, wobei bei Erreichen des gemes- 
senen Ausgangswerts Dv des Systems zum Zeitpunkt tymitto^ty^Tdas Verfahreo unterbrochen wird 
und nach einer diskontinuierlichen Anderung getnaS 



dE / \ dE / \ 



8bx 



dx 



9e 



fortgesetzt wild bis zum nScbsten MeBwert des Systems oder bis zum Erreichen des Anf angszeitpunkts 
(3) durch die Ergebnisse von Schritt (2) werden, bis auf Kcnrekturterme, die gesuchten Gradienten dar- 
gestellt: 



V 







f 








V 




6V(t v ) 



dE 



dE 
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55 



dxo d^O 



(to). 



c) bei dem die Gradienten mimmiert werden und damit das System durch die Abbildung dargestellt winL 
2. Verfahren zur Abbildung eines technischen Systems durch einen Rechner, 

a) bei dem eine Abweichung von gemessenen Ausgangswerten des technischen Systems und von Ausgangs- 
werten der Abbildung dutch folgende Bewertungsfunktion definiert wild: 



e(y(t v ),n v ) = |(y(tv) - n v ) 2 , 



nut 

d 
dt 



x(t) =- f(x(t), u(t), w) und 



y(t) = m(x(t), w) , 

60 

wobei 

t eine Zeitvariable, 
y Ausgabewerte der Modellierung, 
w Systemparameter, 
65 x einen zeitabhangigen Zustandsvektor, 

fQ eine zeitliche Entwicklung der internen Zustande der Abbildung des technischen Systems, 

u Eingabewerte, 

mO eine statische Abbildung, 



6 




V eine Zahlvaria^ 
b) bei dem duich denl 
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gswert des Systems bezeichnen, 
die Beweatungsfimktion dadurch optimiert wild, dafi Gradienten 




dE 
dw 



und 



dE 
dxo 



in folgeoden Scfarinen bestinrmt werden, wobei durch Xq ein Anfangszustandsveklor zum Zeitpunkl \q darge- 
stellt wird: 

(1) das Differentialgleichungssysteni 

A.^t) = f^t), u(t), w ) 



wild von t = tfl bis t = T gelost, unter Vbrwendung des Anfangszustandsvektois Xqt wobei sich ein Endzu- 
standx(T) eigibt, 

(2) ein System von DifferenUalgleichungen 



_d_ dE ^ _ f df(x(t), Ti(t), w)V d£ ( . 
dt dw U ~ ~\ d» J dx 0 V ' 



d dE 



dt dxo 



(t) = \ 



6f (x(t), u(t), w)J dE 



ax 



dXQ 



(t) 



A x(t ) = fftt), u(t), w) 

wobei 0* die transponierte Matrix ist, wild von T rflckwarts bis to gelost, wobei bei Erreichen des genres- 
senen Ausgangswexts Dy des Systems zum Zeitpunkt ^ nut tg < ty S T das Verfahien unleibrochen wird 
und nach einer diskontinuieriichen Andemng gemaB 



dx 0 



dx 0 



an 



dx 



(m(3c(tv> w) - n v ) 



fortgesetzt wird bis zum nachsten MeBwert des Systems oder bis zum Eneichendes Anfangszeitpunkts to, 
(3) durch die Ergebnisse von Schritt (2) werden, bis auf Korreklurterme, die gesuchten Gradienten dar- 
gestellt: 



(m(x(t v > w) - n v ) 



dE 



dE 



dxo dxo 



(to), 



c) bei dem die Gradienten mirrirrriert werden und danrit das System durch die Abbildung daigestellt wird. 
3. Varfahren nach Anspnich 1 oder % bei dem ein Zeitintervall von dem Zeitpunkt to bis zu dem Zeitpunkt T in 
mehrere Teilintervalle aufgeteilt wird, wobei fur jedes Ibilintervall die Verfahrensschritte angewandt werden und 
Einzelergebnisse der jeweiligen Teilintervalle zu einem Gesamiergebnis aufsummiert werden. 
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